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ABSTRACT 

We present observations that were carried out with the Two HUndred Micron Pho- 
tometER (THUMPER) mounted on the James Clerk Maxwell Telescope (JCMT) in 
Hawaii, at a wavelength of 200 /iin (frequency 1.5 THz). The observations utilise a 
small atmospheric window that opens up at this wavelength under very dry conditions 
at high-altitude observing sites. The atmosphere was calibrated using the sky-dipping 
method and a relation was established between the optical depth, r, at 1.5 THz and 
that at 225 GHz: ti.5th z = (95 ± 10) x T225GHz- Mars and Jupiter were mapped from 
the ground at this wavelength for the first time, and the system characteristics mea- 
sured. A noise equivalent flux density (NEFD) of ~65±10 Jy (la Is) was measured 
for the THUMPER-JCMT combination, consistent with predictions based upon our 
laboratory measurements. The main-beam resolution of 14 arcsec was confirmed and 
an extended error-beam detected at roughly two-thirds of the magnitude of the main 
beam. Measurements of the Sun allow us to estimate that the fraction of the power 
in the main beam is ~15%, consistent with predictions based on modelling the dish 
surface accuracy. It is therefore shown that the sky over Mauna Kea is suitable for 
astronomy at this wavelength under the best conditions. However, higher or drier sites 
should have a larger number of useable nights per year. 

Key words: instrumentation: photometers - techniques: photometric - infrared: 
Solar system - submillimetre 



1 INTRODUCTION 

Far-infrared astronomy is largely the preserve of telescopes 
borne on satellites, balloons and aircraft. The Infra-Red As- 
tronomical Satellite, IRAS (Beichman et al., 1988) operated 
at wavelengths as long as 100 fim. The Kuiper Airborne 
Observatory, KAO (Cameron, Bader & Mobley 1971; Har- 
vey 1979), the Infrared Space Observatory, ISO (Kessler et 
al., 1996), and the Spitzer Space Telescope (Werner et al., 
2004) were designed to operate out to 200 /im. A number 
of planned missions such as Herschel and Planck aim to 
cover still longer wavelengths. However, a problem with all 
such airborne and satellite telescopes is that they are un- 



likely to be as large as telescopes that can be built on the 
ground. Consequently, airborne and orbiting telescopes do 
not have the angular resolution of comparable ground-based 
instruments. Therefore, any instruments that can be oper- 
ated from the ground at these wavelengths and mounted on 
large telescopes will automatically have a resolution advan- 
tage. 

Only a limited number of previous ground-based astro- 
nomical observations have been made at far-infrared wave- 
lengths. This is because the atmosphere is largely opaque 
in this regime. However, spectral measurements of the at- 
mosphere, coupled with detailed modelling of atmospheric 
properties, have shown that there are some small windows 
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Figure 1. Plot of normalised crystal response versus frequency 
for an unstressed crystal (lighter line to the right) and a stressed 
crystal (heavier line to the left) in the THUMPER detector ar- 
ray. Note how the peak response has moved to the required 
THUMPER frequency of 1500 GHz. 



that do open up across this range on very high, dry sites 
under good conditions. One such window is that at 200 fim, 
corresponding to a frequency of 1.5 THz. 

Previous ground-based measurements at this frequency 
have only been made spectroscopically. No previous ground- 
based imaging camera has been built to work at this 
frequency. The spectroscopic observations include Fourier 
Transform Spectrometer (FTS) measurements at Pampa la 
Bola in Chile (Matsuo, Sakamoto & Matsushita 1998; Mat- 
sushita et al., 1999). These latter authors derived measure- 
ments relating the optical depths in different atmospheric 
pass-bands, including that at 200 pm. Ground-based FTS 
measurements at Mauna Kea in Hawaii have been carried 
out by a number of groups (e.g. Serabyn & Weisstein 1996; 
Serabyn et al., 1998; Pardo et al., 2004). All have confirmed 
the existence of an atmospheric window at 200 /jm, albeit 
with varying levels of transmission, depending upon the con- 
ditions. 

Paine et al., (2000) carried out similar FTS measure- 
ments at both Mauna Kea in Hawaii and at Chajnantor 
in Chile and measured the atmospheric transmission from 
about 300 to 3000 GHz. They measured the 1.5-THz win- 
dow and found significant transmission in the wave-band 
of interest at both locations. Other work at this frequency 
has included plans for a heterodyne system to operate at 
the South Pole by a group who have made FTS sky mea- 
surements, broadly confirming the above-mentioned results 
(Chamberlin et al., 2003; Gerecht et al., 2003). 

In this paper, we present the results from the first 
commissioning run of the Two HUndred Micron Photome- 
tER (THUMPER) at the James Clerk Maxwell Telescope 
(JCMT). This instrument is a prototype camera built at 
Cardiff University to demonstrate the viability of carrying 
out high-resolution broad-band imaging at 200 /im from the 
ground. 
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Figure 2. Plot of filter transmission versus frequency for the 
complete filter stack in THUMPER. The atmospheric central fre- 
quency corresponds to 1500 GHz. 

2 THE INSTRUMENT 

THUMPER consists of a 7-channel array of photoconduc- 
tor pixel elements arranged on a hexagonal grid (Ward- 
Thompson et al., 2002; Walker et al., 2003; Evans et al., 
2005). This is designed to match the layout of the central 
pixels in the Sub-millimetre Common User Bolometer Array, 
SCUBA (Holland et al., 1999), when placed at the Nasmyth 
focus of the James Clerk Maxwell Telescope (JCMT) on 
Mauna Kea, Hawaii. This is at an altitude of 4092 m, and is 
one of the best sites for sub-millimetre observations in the 
world. 

The THUMPER detectors are photo-conductors made 
from gallium-doped germanium (in the ratio Ge:Ga ~ 10 5 :1) 
crystal semi-conductors, which have been stressed to move 
their peak responsivities into the 200-pm waveband (Walker 
et al., 2003; Walker 2004). Figure 1 shows a typical crystal 
responsivity as a function of frequency. The lighter line to 
the right shows the natural crystal response. The heavier 
line to the left shows the response after the crystal has been 
mounted in the camera and stressed to its optimal value. 
A series of filters matches the band-pass of the instrument 
to the atmospheric window. Figure 2 shows the combined 
band-pass filtering of the system. 

A drawing of one of the stressing blocks is shown in 
Figure 3(a). The crystals are mounted in the integrating 
cavities on the front of the stressing block (on the left of 
Figure 3(a)) and a stressing rod is pushed up into the rear 
of the block (on the right of Figure 3(a)) using a stressing 
screw. This in turn operates a lever which applies the stress 
downwards through the line of crystals by means of further 
stressing rods. 

There are three stressing blocks. The centre block holds 
three crystals and the two side blocks hold two crystals each, 
such that when the three blocks are fixed together, as in 
Figure 3(b), the crystals make a hexagonal arragement. The 
crystal integrating cavities are coupled to the telescope by 
individual feed-horns, as shown in Figure 3(b), with a res- 
olution of 14 arcsec, as for the SCUBA 850-/im array (Hol- 
land et al., 1999). The THUMPER focal plane array can 
also be adjusted to match the SCUBA 450-/im array if de- 
sired, although this mode was not used on this occasion. We 
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Figure 3. The THUMPER focal plane array: (a: top left) a stressing block in which the crystals are mounted; (b: top right) the three 
stressing blocks fixed together to form the focal plane array, with feed-horns on the front of the detectors; (c: lower left) the focal plane 
array mounted on its stand; (d: lower right) the complete array and mount with the filter stack mounted on the front. 



chose not to make the feed-horns diffraction-limited on the 
telescope because the dish surface accuracy is too low to 
provide a significant fraction of the total incident power in 
a diffraction-limited beam at this wavelength. 

The focal plane array is mounted in a stand as shown in 
Figure 3(c) and the filter stack attached to the front of the 
array as seen in Figure 3(d). The whole array is cooled to 
liquid Helium temperature (3.7 K at JCMT) in a long hold- 
time cryostat, built by Queen Mary College Instruments 
(now based at Cardiff University) and Thomas Keating Ltd. 
The system detector quantum efficiency (DQE) was mea- 
sured in the laboratory. Calculations were performed based 
upon the optical set-up of the instrument at the telescope, 
including estimates of the dish surface accuracy. Incorpo- 
rating all of these factors we estimated before going to the 
telescope that the noise equivalent flux density (NEFD) of 
the JCMT-THUMPER arrangement should be -50-70 Jy 
(la Is - Walker 2004). 



3 OBSERVATIONS 

The cryostat was mounted on the right Nasmyth platform 
of the JCMT. Two lenses, made of high-density polyethy- 
lene, were used to bring the beam from the telescope to the 
f/5.75 focus required by THUMPER (Walker et al. 2003). 



The instrument was commissioned during Director's Discre- 
tionary Time in two 4-night observing runs on 2005 March 
20-23, and 2005 April 6-9. Optical alignment and instru- 
ment setup were carried out prior to this. Data acquisition 
was carried out via a stand-alone computer operating the 
Lab VIEW software system (Laboratory Virtual Instrument 
Engineering Workbench) . Only one of the scheduled nights 
was good enough for astronomical observations at this wave- 
length, 2005 April 9, when our observing was scheduled from 
01.30 to 09.30 Hawaiian Standard Time (UT 11.30 to 19.30). 
Data from that night are presented in this paper. 



4 THE 200-MICRON SKY 

The 200- /jm (1.5 THz) atmospheric window is much nar- 
rower in wavelength than the more commonly used sub- 
millimetre windows, and has much lower transmission. We 
modelled this window, along with neighbouring windows, us- 
ing the Cardiff-developed software package FETCH (Araujo 
et al., 2001; Hayton et al., 2005). This is a fast, highly flexi- 
ble line-by-line, layer- by- layer, radiative transfer model that 
uses the latest High-resolution Transmission (HITRAN) 
molecular spectroscopic database (2004 version) to calcu- 
late the atmospheric transmission as a function of frequency 
(Rothman et al., 2003; 2005). We found that the window 
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Figure 4. Plot of atmospheric transmission versus frequency as 
we have modelled, for 0.5mm PWV. Note the more frequently 
used sub- millimetre window at 850 GHz (350 (im). The 1.5-THz 
window is beginning to open up in these conditions. 
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Figure 5. Plot of a typical set of sky-dip data. Normalised emis- 
sion is plotted against zenith distance. The data-points arc shown, 
along with estimated error-bars, joined by a dashed line. The solid 
line is a fit to the data used to obtain T\,^tHz ■ See text for details. 

extends from 1428 to 1588 GHz (189-210 /im), with peak 
transmission at 1486 GHz (202 /im). The reason why the 
peak wavelength is not centred in the window is that there 
is an absorbtion line directly in the centre of the window. 

Figure 4 shows the result of the detailed model calcula- 
tions of the transmission across the whole of this part of the 
spectrum, assuming that the atmospheric precipitable water 
vapour (PWV) content is 0.5 mm. The more familiar window 
at 850 GHz (350 pm) can be seen, along with the small trans- 
mission window just opening at 1.5 THz (200 /im). The peak 
transmission is low, but nonetheless useable for part of the 
time. This calculation is consistent with the measurements 
of Paine et al. (2000), under better atmospheric conditions 
in Chajnantor, Chile - see their figure 3. 

On the occasions when the sky was totally opaque we 
measured the relative responsivity of the seven channels by 
measuring the sky emission at zenith. This is essentially a 
flat, extended source, and allowed us to make a flat-field of 
our focal plane array. This remained remarkably constant, 
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Figure 6. Plot of ti.qtHz versus TiihGHz- The data arc shown 
as crosses. The solid line is a least-squares fit to the data. The 
best fit straight line that we find is ti.sthz = (95 ±10) X T225GHz- 

and was also consistent with our sensitivity measurements 
made in the laboratory. 

During the times when we could carry out astronom- 
ical observations we calibrated the sky transmission using 
the method of sky-dipping. This entails measuring the sky 
emission at several positions between zenith and horizon, 
and modelling the profile as a function of zenith distance, z. 
The emission is predicted to have the form (1 — e~ Tsecz ), and 
hence a fitting routine can be used to calculate the optical 
depth, r. 

Figure 5 shows the result of a typical sky-dip. On the 
x-axis is plotted zenith distance z from 90° to 0°, while the 
y-axis gives the emission measured from the sky at a given 
value of z, normalised to the peak emission at 90° . The data- 
points with error-bars are shown connected by a dashed line, 
while the solid line shows the model fit to the data for an 
optimised value of the optical depth, r. 

Throughout our observations, the JCMT water vapour 
monitor (WVM) measured the atmospheric opacity at 
183 GHz and from this calculated the opacity at 225 GHz 
(for historical reasons) in the standard way in which it 
records these data (e.g. Archibald et al., 2002). The WVM 
operates via a pick-off mirror at the edge of the JCMT field 
of view. It updates its estimate of the 225-GHz opacity ev- 
ery 1.2 seconds. Consequently, at the time of each sky-dip 
we know the value of the 225-GHz opacity, so we can plot it 
against our measured 1.5- THz opacity. 

Figure 6 shows the resultant plot of ti.sthz versus 
T225GHZ ■ The data are shown as crosses. Our sky model cal- 
culations (see above) predict a linear relation between these 
two frequencies over the observed range of values. The solid 
line is the resulting least-squares fit to the data. The best 
fit straight line that we find is: 

Tl.STHz = (95 ± 10) X T225GHZ- 

This result is consistent with our model predictions and can 
also be compared with those obtained by previous workers. 
For example, Matsushita et al. (1999), found a relation of 
the form ti.bthz = (105 ± 32) x T225GHZ for the observing 
site at Pampa la Bola in Chile. 

These latter authors used an FTS to measure the at- 
mospheric opacity across the different wavebands, whereas 
we are comparing our data to WVM data. Nevertheless, re- 
cent work shows that provided all relevant effects are taken 
into consideration, the results from FTS and WVM mea- 
surements generally agree well (Pardo et al. 2004). Hence 
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Jupiter at 200 microns 



Mars at 200 microns 
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Figure 7. First-light astronomical image taken with THUMPER: 
A 200-^iin isophotal contour map of Jupiter. The field of this 
image roughly corresponds to the size of the THUMPER array. 
The centre of the array is at (0,0). Contour levels arc at 50, 70 
& 90% of peak. The fact that the image appears slightly non- 
circular shows that part of the error-beam of the telescope may 
be affecting the data. The optical alignment is also slightly offset. 



R.A. offset (arcsec) 

Figure 8. A 200-^im isophotal contour map of Mars taken by 
THUMPER on JCMT. The contour levels are at 50, 60, 70, 80 
& 90% of peak. Once again the array centre is at (0,0). The 
alignment is improved in R.A. The extension to the north is the 
telescope error beam, which probably also extends to the south, 
although this was not mapped. 



our result is fully consistent with previous work. Therefore 
we can use this relation to calibrate our subsequent mea- 
surements. 



5 PLANET MAPS 
5.1 Jupiter 

The first astronomical object we imaged was Jupiter. Fig- 
ure 7 shows our 200-/im map of Jupiter. It was constructed 
from two consecutive maps of Jupiter, taken immediately 
one after the other, over the airmass range 1.25 to 1.27. 
The two maps both show the same structure, and so they 
were co-added to increase the signal-to-noise ratio. We note 
that the source is not exactly centred in the image, being 
misaligned by ~30 arcsec. The telescope absolute pointing 
accuracy (based on measurements made with other instru- 
ments) at this time was <3 arcsec. This implies that our 
system alignment was not perfect. This is not surprising 
given that this was the first time that THUMPER had been 
mounted on JCMT and this was the first-light astronomical 
image. 

The flux density of Jupiter at this wavelength can be 
predicted from planetary modelling (Griffin et al., 1986; Or- 
ton et al., 1986; Griffin & Orton 1993). Based on this, we 
expect the peak flux density of Jupiter to be 34 kjy/beam. 
During the observations the measured value of T225GHZ 
was 0.0625. Using our relation above, between 1.5THz and 
225GHz, this corresponds to a value of ti.sthz of ~5.95. 

Therefore the observed flux density at the telescope 
is predicted to be 21 Jy/beam. The total integration time 
per point of the two co-added maps was 100 seconds. The 
peak was detected at a level of ~3.3 a. Therefore we cal- 
culate from the Jupiter data that the noise equivalent flux 
density (NEFD) of the JCMT-THUMPER combination is 
~63±10 Jy (la Is). 

We note that our detection is not at the 5-cr level. How- 



ever, we believe it is a real detection for a number of rea- 
sons: the source was seen in several pixels simultaneously; 
the structure in the two maps that we took was the same; 
the source appeared in the same place in both maps; and 
the NEFD we calculate from the measurements is consis- 
tent with that predicted from laboratory measurements of 
the detector system. In re-gridding the map onto an RA- 
Dec grid as shown in Figure 7 some smoothing naturally 
occurred, as a pixel scale of 5 arcsec was adopted. 

The full- width at half-maximum (FWHM) of Jupiter in 
our map is ~46x38 (±8) arcsec. At the time of the obser- 
vations, the diameter of Jupiter was 42 arcsec. When con- 
volved with our 14-arcsec beam, this becomes 44.2 arcsec. 
Hence our observations are consistent with this. However, 
the non-circularity visible in the map shows that we are also 
detecting part of the telescope error beam, which is pre- 
dicted to be significant at this wavelength. We could not 
trace the error beam further in these data due to their low 
signal-to-noise ratio. 

It is also possible that if we did not have the instrument 
exactly in focus then this would also contribute to the appar- 
ent non-circularity of the beam. We focussed the instrument 
using our model of the telescope and our calculation of the 
optimal focal position. However, we did not have time to 
check the focus on-source before Jupiter began to set. 

5.2 Mars 

Later in the night we imaged Mars as it was rising, from 
airmass 1.30 to 1.29. Figure 8 shows our map of Mars. Once 
again the source is not centred, although we have managed 
to improve the alignment in R.A. at least. The flux density 
of Mars at this wavelength can be predicted from plane- 
tary modelling in the same way as that of Jupiter quoted 
above (Griffin et al., 1986; Orton et al., 1986; Griffin & Or- 
ton 1993). Based on this, we expect the peak flux density 
of Mars to be 8.1 kjy/beam. During the observations the 
measured value of T^ihGHz was 0.0455. Using our relation 
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above, between 1.5 THz and 225 GHz, this corresponds to 
a value of ti.sthz of ~4.32. 

Therefore the observed flux density at the telescope is 
estimated to be 30 Jy/beam. The integration time per point 
of the map was 50 seconds. The peak was detected at a 
level of ~3.2 a. Therefore we calculate from the Mars data 
that the noise equivalent flux density (NEFD) of the JCMT- 
THUMPER combination is ~66±10 Jy (Iff Is). 

We note once again that our detection is not at the 5- 
<7 level. However, we believe it is also a real detection for 
similar reasons to those quoted above: the source was seen 
in several pixels simultaneously; and the NEFD we calcu- 
late from the measurements is consistent with that predicted 
from laboratory measurements of the detector system and 
with that seen in the Jupiter data. 

The map of Mars shows a slightly different morphol- 
ogy from that of Jupiter. The map shows some evidence 
that the brightest region is more centrally peaked, and that 
this central peak sits on an extended plateau. Conversely, 
Jupiter is better fitted by a single gaussian. For Mars the 
centrally peaked core is of order ~15 arcsec across, with a 
more extended lobe to the north. The core is consistent with 
the FWHM that would be expected from Mars, since at the 
time of the observation the diameter of Mars was 6 arcsec. 
When convolved with our beam this produces a gaussian of 
FWHM 15.2 arcsec, consistent with the image in Figure 8. 

The extended emission to the north is most likely to be 
the error-beam of the telescope. The level of the extended 
error-lobe is roughly two-thirds of the magnitude of the main 
beam. This may also extend to the south and south-east, but 
our map does not extend far enough in these directions to 
say. Once again the instrument focus may be adding to the 
problem. By the time we repeated the map the shift was 
coming to an end and the sun had risen. Consequently, the 
conditions worsened (the sky noise level increased) and we 
did not detect Mars again in day-time. Our mean NEFD 
measured on the planets is therefore ~65±10 Jy (la Is). 



6 THE SUN 

Once the Sun had risen we pointed the telescope at the 
Sun, when it was in the airmass range 1.44 to 1.43, and 
detected it clearly in all channels, at high levels of signal- 
to-noise ratio of up to 560 a. Obviously, we did not map 
the full extent of the Sun, but rather used it as a bright, 
uniform, extended black-body source, and chopped across 
the limb of the Sun. We estimate that the emission from 
the Sun is ~1.1 MJy/beam. It was at ~1.4 airmasses during 
our observation, and T225GHZ was 0.047. We used 5-second 
integrations, so we estimate our NEFD on the Sun to be 
~9 Jy (Iff Is). 

This is a factor of ~7 better than our NEFD estimated 
from the essentially point-like planet Mars. This implies that 
only roughly one-seventh (~15%) of the total power incident 
upon the JCMT dish at 200 /jm is focussed into a 14-arcsec 
central beam. We had predicted that, due to the surface 
inaccuracies of the JCMT dish (which was not designed to 
operate at this high a frequency), there would be significant 
power in the side-lobes. 

The measurement of the surface accuracy of the dish 
which is closest in time to our observations was taken on 



2005 February 13, at which time the dish surface accuracy 
was found to be 23.8 /im (Wouterloot 2005). Using the stan- 
dard calculation of Ruze efficiency we estimate that a dish 
of this accuracy would concentrate ~11% of the total power 
into a 14-arcsec beam. This is consistent with the ~15% we 
estimate here based on the ratio of NEFDs measured on 
Mars and the Sun. Hence we see that all of our observations 
are self-consistent. 



7 CONCLUSIONS 

We have successfully commissioned the THUMPER camera 
at JCMT. We have demonstrated that 200-/im astronomy is 
possible from the ground. We have taken the first ground- 
based images of Jupiter and Mars at this wavelength. We 
have calibrated the 200-/im sky, and found a relation be- 
tween opacities of ti.sthz = (95 ± 10) x T225GHz- This is 
consistent with previous measurements, and also with our 
modelling of atmospheric transmission across this section of 
the electro-magnetic spectrum. We estimate the NEFD of 
the JCMT-THUMPER system is ~65 Jy (Iff Is). We find 
that the error beam of JCMT is substantial at this wave- 
length, as predicted, with only ~15% of the power appear- 
ing in the central 14 arcsec. We note that a telescope with 
a greater surface accuracy would improve upon this num- 
ber, and that other observing sites in even higher or drier 
locations may provide a larger number of useable observing 
nights per year. Based on its current NEFD, it seems that 
THUMPER could be used to study bright, high-mass, star- 
forming regions at higher angular resolution than has been 
previously possible at this wavelength. 
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